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ABSTRACT 

In order that magnetic flux be confined within the solar interior for times comparable to 
the solar cycle period it has been suggested that the bulk of the solar toroidal field is stored 
in the convectively stable overshoot region situated beneath the convection zone proper. Such a 
magnetic field though is still buoyant and is therefore subject to Rayleigh-Taylor type instabilities. 
In this paper we consider the model problem of an isolated region of magnetic field embedded in a 
convectively stable atmosphere. The fully nonlinear evolution of the two-dimensional interchange 
modes is studied, thereby shedding some light on one of the processes responsible for the escape 
of flux from the solar interior. 

1. INTRODUCTION 

Despite considerable strides in recent years in both dynamo theory and the study of magne- 
toconvection, several problems remain regarding certain key features of the sun’s magnetic cycle. 
One of these concerns the location of the sun’s primary toroidal field and the means by which 
this can escape from the interior, eventually to appear at the surface to form sunspots or smaller 
bipolar active regions. The inherent buoyancy of isolated magnetic flux tubes, together with the 
unstable convective motions, will carry tubes of equipartition strength through the entire convec- 
tion zone on a timescale of about a month (see, for example, Parker, 1975; Moreno-Insertis, 1983). 
This presents certain difficulties for keeping the requisite amount of magnetic flux within the con- 
vection zone for timescales comparable to that of the solar cycle. To overcome this difficulty it 
has been suggested (Spiegel and Weiss, 1980) that the bulk of the toroidal field is contained in 
the convectively stable overshoot region situated beneath the convection zone proper. Of course 
such a field alters the density stratification, causing the atmosphere to become top-heavy — in 
consequence, Rayleigh-Taylor instabilities can ensue. 

In this paper we consider the instabilities which can result when an isolated region of 
magnetic field is embedded in a convectively stable atmosphere. It should be noted that this is 
a very different situation to that of an isolated flux tube in a field-free gas. In particular, our 
initial static state possesses a continuous temperature distribution and is a bona fide solution 
of the hydrostatic equation (see Section 2) — for an isolated flux tube in thermal equilibrium 
with its surroundings, on the other hand, there is no static state and the tube is forced to rise. 
Our calculation is both fully compressible and nonlinear but, for simplicity, we shall here only 
consider the two-dimensional interchange modes (no bending of the field lines) and shall ignore 
the effects of rotation. The account given is, of necessity, rather brief — a much more detailed 
version, which also treats rotational effects, is currently in preparation. 

2. MATHEMATICAL FORMULATION 

We consider a plane parallel layer of conducting fluid of thickness d. Within the layer 
the magnetic field B and gravity g are constrained to be mutually orthogonal for all times. The 
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fluid flow is in the xz plane only, perpendicular to the magnetic field. The computational domain 
can be thought of as a section of a meridional cut near the equator. The field is then toroidal 
and the flows are axisymmetric. To simplify the problem the thermal conductivity K, the shear 
viscosity /z and the magnetic diffusivity rj are assumed to be constants; furthermore the plasma 
obeys the perfect gas law with constant specific heats. By choosing the layer depth d and the 
sound crossing time dj y/p/p as the units of length and time respectively the evolution equations 


can be expressed in dimensionless form as follows: 

P = PT, ( 1 ) 

d t p + V • pu = 0, (2) 

d t B + V • {Bu) = rC k V 2 B, (3) 

dtpu + V • puu = -V(p+ ^-) + 0{m + l)pz + C fc a(V 2 u + ^V(V -u)), (4) 


p(d t T + u • VT) + (7 - l)pV • u = 7 C fc V 2 T + (7 - 1 )C k {ad i u j <f>ij + r(V£) 2 //?), (5) 


where <f>ij = d{Uj + djm — |<5{y V • u. 

The six dimensionless parameters which appear explicitly are defined by: 
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where A is the initial temperature gradient. 

The formulation of the problem is completed by imposing boundary conditions at the 
top and bottom and by requiring that all variables be periodic in the horizontal direction. The 
boundary conditions at z = 0 (top) and z = 1 (bottom) are chosen so as to minimise their effects 
in the interior; they are: 


T = 1 at z = 0 and T = 1 + 0 at z = 1, 


pw = d z u = d z B = 0 at z — 0,1. 

If the magnetic field has a top hat profile then the equations admit a piecewise-polytropic 
static solution. Small random disturbances (typically about 1 %) of this equilibrium are used 
as initial conditions for the fully nonlinear calculations. The stability of the static solution to 
infinitesimal disturbances is determined by linear theory. 


3. THE LINEAR REGIME 

The linearised versions of Eqs. (1) to (5) are separable in x and t with solutions of the form 
f(x,z,t ) = f(z)e 8t coslx (or sin£x). We are then left with a set of coupled ordinary differential 
equations (involving only z) with the growth rate s as the eigenvalue. For small values of r 
instability always occurs via an exchange of stabilities. We find that there is always one dominant 
unstable mode and that this has grown appreciably before any of the harmonics lose stability. 
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As expected, increasing the strength of the magnetic field (reducing /3 ) leads to the instability 
becoming more vigorous (see figure l). The figure also shows that as the field is increased then 
the mode of maximum growth rate moves to larger wavenumbers (smaller cell widths). In the 
absence of all diffusion the growth rate is maximised for infinitesimally thin cells (£ — ♦ oo); 
however, dissipation acts most effectively on such cells and therefore a balance is struck at a finite 
cell size. For small /? the instability is sufficiently strong that fairly small cells can overcome the 
dissipative effects and grow most rapidly. This behaviour is analogous to that which occurs in 
the standard Rayleigh-Taylor instability of two superposed viscous fluids (Chandrasekhar, 1961, 
p.446). A more surprising feature is that the dependence of the mode of maximum growth rate 
on the depth of the magnetic layer is very weak (see figure 2). This would appear to be because 
the instability occurs primarily at the upper interface of the magnetic field and that the depth of 
the field below that is not of fundamental importance. 



Figure 1 Contours of the growth rate for 
/? versus i. o = .1, r = .05, Ck = .05, 
m = 1 . 6 , 0 = .1 



Figure 2 Contours of the growth rate for h, 
the depth of the magnetic region, versus 
L a = t = .1, Ck = -05, m = 1.6, 0 = 1 


4. THE NONLINEAR REGIME 

The evolution of a uniform magnetic field confined initially to the region 0.6 < z < 0.8 is 
shown in figure 3. The instability due to magnetic buoyancy leads to an initial break-up of the 
layer and, eventually, to the spreading of magnetic field over the entire domain in a time (a few 
Alfven crossing times) short compared to the diffusion time. 

In the initial phases of the instability the magnetic layer’s upper interface forms the familiar 
mushroom-shaped structures associated with Rayleigh-Taylor instabilities (see, for instance, Daly, 
1967). In all cases considered the density difference between the magnetic and field-free regions 
was such that mushrooms were formed rather than bubbles and spikes. Associated with the wings 
of the mushrooms are regions of intense vorticity which, once most of the available potential energy 
has been released, play a key role in the subsequent evolution of the layer. The characteristic 
wrap-around structures visible in figures 3(c) to 3(f) are typical of flows generated by random 
distributions of vortex pairs. Some of the vortex-vortex interactions are strong enough to do, 
albeit locally, negative buoyancy work and to produce the downward displacement of regions of 
high magnetic field even though these regions are lighter than their surroundings and would have, 
in the absence of any other flow, risen upwards. 
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Figure 3 Time evolution of an unstable magnetic layer, o = .1, r = .01, C* — .01, m — 1.6, 

6 = 2, /? = 1. The simulation times for the frames are respectively: (a) 1.9, (6) 3.5, (c) 5.1, 
(d) 6.7, (e) 8.3 and (/) 13.1. 
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